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Abstract 



The space-time evolution of the hot and dense matter formed after the 
colhsions of heavy nuclei at ultra-relativistic energies is investigated using 
^ . (3+1) dimensional hydrodynamical models. The effects of the spectral shift 

! of the hadronic properties are incorporated in the equation of state (EOS) of 

' the evolving matter. In- medium shift of hadronic properties are considered for 

Quantum Hadrodynamics (QHD) and universal scaling scenarios. It is found 
! that the EOS for the hadronic matter for universal scaling of hadronic masses 

(except pseudoscalar) is similar to the recent lattice results. We observe that 
the space-time volume of the hadronic matter at the freeze-out is considerably 
different from the one when medium effects on the hadrons are ignored. The 
sensitivity of the results on the initial radial velocity profile is investigated. 
The transverse mass spectra of pious and protons of NA49 collaboration are 



P3 . analyzed. 



^ : 1 Introduction 



One of the main motivations to study the nucleus-nucleus collisions at ultra-relativistic 
energies is to create a very hot and dense system of strongly interacting matter, a 
situation conducive for the formation of Quark Gluon Plasma (QGP) [1]. It is not 
easy, however, to obtain information about the early stage of this matter because 
of the very short lifetime of the dense system. Electromagnetically interacting par- 
ticles (photons and lepton pairs) are considered to be the ideal probes for the early 
state of the matter. On the other hand, observations through the transverse mass 
(mr) spectrum of the hadrons provide valuable informations about the situation 
when the thermal system freezes out [2, 3, 4, 5, 6], i.e. the stage when the system 
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disassembles to individual hadrons. Thus the study of the hadronic spectra gives 
snapshots of the later stages of the nuclear collision dynamics. 

It has been emphasized recently that the properties of the hadrons will be mod- 
ified due to its interactions with the particles in the thermal bath. It is also worth 
emphasizing here that as yet it has not been possible to explain the enhanced dilep- 
ton production measured by the CERES/NA45 collaborations [7] in the low in- 
variant mass region (below the p-peak) without the in-medium modifications of the 
vector mesons [8]. We have recently analyzed [9] the WA98 photon data [10] by 
incorporating the in-medium modifications of hadrons. In the light of these observa- 
tions we study the hadronic spectra in a scenario where the hydrodynamic evolution 
contains the in- medium modifications of hadrons through the parametrization of the 
equation of state (EOS). In particular we will show here that the NA49 [11] hadronic 
spectra can be explained by the same initial and freeze-out conditions with the EOS 
which explain the photon and dilepton data at CERN SPS energies. We consider 
two possible scenarios: (i) nucleus -|- nucleus — > QGP — > hadrons and (ii) nucleus 
+ nucleus — ^ excited hadronic matter — > hadrons (but the properties of hadrons 
in the thermal bath are different from the vacuum due to its interaction with the 
particles in the thermal bath), to investigate the evolution of the matter formed 
after the collisions. 

In the next section we review the hydrodynamical model and the EOS used in 
the prersent work. In section 3 we present the results of our calculations and section 
4 is devoted to summary and discussions. 

2 The Single-particle Spectra in the Hydrody- 
namical Model 

We will assume that the produced matter reaches a state of local thermodynamic 
equilibrium after a proper time, ~ 1 fm/c [12]. The system continues to expand 
in space and time till the freeze-out undergoing a phase transition to the hadronic 
matter in the process, if QGP is formed initially. The conversion of thermal energy 
of the system into the collective flow continues as long as coUisions are sufficiently 
frequent. Eventually the density of the particles decreases and correspondingly their 
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mean free path increases. When the mean free path becomes of the order of the size 
of the system, hadrons act as non-interacting free particles. The detailed description 
of this decoupling transition, or freeze-out, is complicated. A simple algorithm is 
obtained if one assumes that the break-up occurs when the temperature drops down 
to a given decoupling temperature Tp ; this temperature, which is typically of the 
order of the pion mass, defines a space-time surface T(r, t) — Tp, usually referred to 
as the freeze-out surface (FS), on which the system ceases to behave collectively. A 
fluid element which crosses this surface liberates particles which, in the rest frame 
of the fluid element, have a thermal distribution at a temperature Tp. Knowing the 
velocity of the fluid element on the decoupling surface, one can then determine the 
momentum distribution of the final particles. The observed momentum distribution 
will thus be characterized by the temperature and collective velocity of the system at 
the FS. For given initial conditions and EOS we solve the hydro dynamical equations 
to obtain the FS, determined by the condition, T(r, r) = Tp. 

The expansion of the system is described by the energy-momentum conservation 
of an ideal fluid : 

where e is the energy density, P is the pressure measured in the frame co-moving 
with the fluid, and is the fluid four-velocity. 

With the assumption that the system undergoes a boost-invariant longitudinal 
expansion along the z axis [12] and a cylindrically symmetric transverse expansion 
the hydrodynamic Eqs. 1 reduce to 

S^yoo ^ la^(^7^oi) ^ 1 ^yoo + p) ^ (2) 
r T 
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where u^^ — "y{t/T, Vr, z/t) 

Eqs. (2) and (3) are solved by the relativistic version of the flux corrected trans- 
port algorithm [13] with the following initial energy density profile, 

'^'■^'''^^ l + exp[(r-i?^)/5]- 
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where Ra is the nuclear radius and 6 is the surface thickness parameter. We have 
taken 5 = 1/2 in our calculations. We will study the sensitivity of the flow with the 
following two kinds of initial velocity profile [14, 15], 



T 

Vr{ri,r) = vo— (5) 



and [13] 



Vr{Ti, r)=vo 1 - — ^ (6) 

[ l + ex.p[{r - Ra)/S]\ 

The FS, which is required as an input to evaluate the single particle spectra is 

obtained by solving the Eqs. (2) and (3) with the initial conditions (4) and (5) (or 

eq.6) for a given EOS. 

The single particle transverse momentum distribution in the hydrodynamical 

model is given by the well-known Cooper- Frye formula [16], 

where f{r,p) is the Boltzmann distribution (quantum statistical effects are neglected 
here) 

f{r,p) = exp[-(p.'u - i^)/T]. (8) 

T, II and u'^ are the temperature, chemical potential and four-velocity of the fire- 
cylinder respectively. The three dimensional space-time surface a, whose surface 
elements are specified in terms of the four vector da^ — {d^r, dtdS), is determined 
by the freeze-out condition, T(r, r) = Tp as mentioned before. For a system with 
cylindrical symmetry and boost invariance along the longitudinal direction the single 
particle momentum (p^) distribution is given by [14, 17, 18], 
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where Tpir) refers to the freeze-out time which, in general, depends on r and yt is 
the transverse rapidity. 
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2.1 Equation of state 

The set of hydrodynamic eqs. (2) and (3) are not closed by itself; the number of un- 
known variables exceeds the number of equations by one. Thus a functional relation 
between any two variables is required so that the system becomes deterministic. The 
most natural course is to look for such a relation between the pressure P and the 
energy density e. Under the assumption of local thermal equilibrium, this functional 
relation between P and e is the EOS. Obviously, different EOS's will govern the hy- 
drodynamic flow quite differently and as far as the search for QGP is concerned, 
the goal is to look for distinctions in the observables due to the different EOS's 
(corresponding to the novel state of QGP vis-a-vis that for the hadronic matter). 
It is thus imperative to understand in what respects the two EOS's differ and how 
they affect the evolution in space and time. 

A physically intuitive way of understanding the role of the EOS in governing 
the hydrodynamic flow lies in the fact that the velocity of sound = {dP/de)s 
sets an intrinsic scale in the hydrodynamic evolution. One can thus write a simple 
parametric form for the EOS: P — c^(T')e. Clearly, inclusions of interactions may 
drastically alter the value of c^. In the present work MIT bag model equation of 
state is assumed for the QGP where the energy density and pressure are given by 



The effective degrees of freedom in QGP, qq — 37 for two flavours. The entropy 
density sq is given by Sq = 2gQ{7i^/45)T^. 

As the hydrodynamic expansion starts, the QGP begins to cool until the tem- 
perature drops down to the critical temperature Tc. At this instant, the phase 
transition to the hadronic matter starts. Assuming that the phase transition is a 
first-order one, the released latent heat maintains the temperature of the system at 
the critical temperature T^, even though the system continues to expand. The cool- 
ing due to expansion is compensated by the latent heat liberated during the process. 
We neglect the scenarios of supercooling or superheating and any other plausible 
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explosive events. This process continues until all the matter has converted to the 
hadronic phase with the temperature remaining constant at T = Tg. From then on, 
the system continues to expand, governed by the EOS of the hot hadronic matter 
till the freeze-out temperature Tp at the proper time tf- Thus the appearance of 
the so called mixed phase at T = Tc, when QGP and hadronic matter co-exist, 
is a direct consequence of the first order phase transition. Apart from the role in 
QGP diagnostics, the possibility of the mixed phase affects the bulk features of the 
evolution process also. 

In the hadronic phase one must consider the effect of the presence of heavier 
particles and the change in their masses due to finite temperature effects. The 
ideal limit of treating the hot hadronic matter as a gas of pions originated from 
the expectation that in the framework of local thermalization the system would 
be dominated by the lowest mass hadrons while the higher mass resonances would 
be Boltzmann suppressed. Indirect justification of this assumption comes from the 
experimental observation in high energy collisions that most of the secondaries are 
pions. Nevertheless, the temperature of the system is higher than during a major 
part of the evolution and at these temperatures the suppression of the higher mass 
resonances may not be complete. It may therefore be more realistic to include higher 
mass resonances in the hadronic sector, their relative abundances being governed 
by the condition of (assumed) thermodynamic equilibrium. We assume that the 
hadronic phase consists of tt, p, a;, 77, oi mesons and nucleons. The nucleons and 
heavier mesons are expected to play an important role in the EOS particularly, in 
a scenario where mass of the hadrons decreases with temperature. 

The energy density and pressure for such a system of hadrons are given by 

d^pENfFD{EN,T) (12) 

2 

d'p^fFDiE^.T) (13) 

where the sum is over all the mesons under consideration and N stands for nu- 
cleons. The in-medium mass of the hadrons {H) enters through the relation. 
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Eh = yp"^ + iTT'H (the asterisk denotes the effective mass in the medium). The 
entropy density is then given by 

eH + Pn 



4aeff(T)T3 = 4^5,ff(m*(T),T)T3 



sh 



T 



(14) 



where Qes is the effective statistical degeneracy. Thus, we can visualize the finite 
mass of the hadrons having an effective degeneracy ges{m*{T),T). We consider the 
effects of in-medium mass on the EOS both for the universal scaling scenario [19] 
and the Quantum Hadrodynamic (QHD) scenario. According to the universal scal- 
ing scenario the hadronic masses (except pseudoscalars) approaching zero near the 
critical temperature as follows [19], 



2\\ 



mymH = (1 - T'iTt) 



(15) 



Results for various values of the exponent A will be discussed below. For a detailed 
discussions on the in-medium masses of hadrons in QHD interaction we refer to 
[20, 21, 22]. The temperature dependence of the mass of nucleon (m^), p [ni*^ and 
a; (m* ) in QHD are parametrized as follows: 
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The effective mass of the ai meson, m*^ has been estimated from m* by using 
Weinberg's sum rule [23]. 



3 Results 

In Fig. 1 the variation of the effective degeneracy (the co-efficient of e/T^ differs from 
the effective degeneracy by a multiplicative factor of 7r^/30) is plotted as a function 
of temperature. The rapid increase of the effective degeneracy near the critical 
temperature as per lattice QCD calculations [24] is reasonably well described when 
the effective masses vary according to the eq.(15) (see also [25]). 
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Figure 1: The energy density e in unit of for various equations of state used 
in the present work are plotted as a function of temperature in unit of the critical 
temperature, Tc. The filled circles denote the lattice results [24]. The dashed line 
indicates results with vacuum masses. 

We solve the following equation to estimate the initial temperature, 

dN^ 45C(3) 2. /^^^3 no^ 
= ^^7ri?^4aeif(Ti)7;Vi (19) 

where aefr(Tj) = {n"^ / 90) gcs{m*{Ti),Ti). The change in the expansion dynamics 
as well as the value of the initial temperature due to medium effects enters the 
calculation through the effective statistical degeneracy. 





QGP 


vacuum 

mass 


QHD 


Scaling 
Tc=200 MeV 


dN/dy 


Ti (MeV) 


Ti (MeV) 


Ti (MeV) 


Ti (MeV) 


700 


196 


245 


220 


205 



Table I: Initial temperatures for QGP and hadronic initial states. 

We consider Pb + Pb collisions at CERN SPS energies. Taking dNyr/dy = 700 
for Pb + Pb coUisions and Ti — 1 fm/c, Eqs. (14) and (19) are solved self consistently 
to get the initial temperatures for different mass variation scenarios shown in Table 
I. For a given dN/dy and Tj the value of Ti is higher in case of vacuum masses 
(smaller Og//) compare to the universal scaling scenario (higher Og//). 

Next we consider the constant energy density contours at the FS. The value of 
Tp is taken as 120 MeV in the present calculations. Eqs. (4) and (5) are used for the 
initial energy density and velocity profiles respectively. In case of QGP formation 
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Energy Density Contours at the freeze-out surface 
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Figure 2: Constant energy density contours at the FS (in r — r plane) for initial 
velocity profile of Eq. 5 with Vq — 0. Left (Right) : QGP (hadronic) initial state. 

the initial energy density is 2.6 GeV/fm^ corresponding to an initial temperature 
of 196 MeV. In this case the in-medium effects on the hadronic phase formed after 
the phase transition is taken into account. The energy density at the freeze-out 
point, ep (oc aeff(Tp)Tp) is different for different hadronic model, ep is larger for 
the universal scaling scenario with A = 1/2 compared to the case of the hadronic 
gas with vacuum masses, because for a given Tp, a^ff is larger in the former case. 
As a consequence of this the FS is smaller in universal scaling scenario compared 
to the vacuum scenario. This is demonstrated in the left panel of fig. 2. The 
FS corresponding to QHD and the vacuum are indistinguishable because of the 
small medium effects in QHD especially at the late stage of the evolution, where 
the temperature is lower. In the right panel of fig. 2 we display the constant ep 
surfaces for hadronic initial state. The velocity of sound in case of universal scaling 
is smaller than the scenario when medium effects are ignored. As a result the FS for 
the universal scaling scenario is larger (for a given Tp) as compared to the vacuum 
mass scenario, because of the slower cooling of the system in the former case. 

In fig. 3 the results for non-zero initial radial velocity on the surface of the 
cylinder is depicted. Because of the rapid expansion of the system the freeze-out 
surface is smaller {i.e. the life time is shorter) in this case as compared to the case 
where the initial radial velocity is zero at the surface. However, the qualitative 
shape of the curves are same as before. 

Fig. 4 shows the FS when a different initial velocity profile, Eq. 6 is used with 
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Figure 3: Same as fig. 2 for vq = 0.2. 




Figure 4: Same as fig. 2 with initial velocity profile of Eq. 6 for vq — 0.2. 
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Constant Velocity Contours 



Constant Velocity Contours 




r (fm) r (fm) 

Figure 5: The constant velocity contours for the initial velocity profile of Eq. 5 with 
vq — 0. Left (Right): QGP (hadronic initial state for vacuum masses of hadrons). 

Vq — 0.2. It is clear that here the expansion is slower than the case in fig. 3, 
indicating that the radial velocity proportional to radial co-ordinate gives rise to 
stronger fiow. 

In fig. 5 we show the constant velocity contours for QGP and hadronic initial 
state when the medium efi^ects on the hadrons are ignored. The initial velocity 
profile is given by Eq. 5 with vq = 0. Introduction of the medium effects does not 
change the qualitative behaviour of the contours as shown in fig. 6. In case of the 
QGP formation scenario the introduction of the medium effects has negligible effects 
on the space time evolution of the system (left panel of figs. 5 and 6). However, 
in case of the hadronic scenario (right panel of figs. 5 and 6) a visible quantitative 
change is evident. In fig. 7 the effects of non-zero initial radial velocity are shown, 
a quantitative difference with the results of the previous figure indicates a large 
transverse flow in this case. 

Now we discuss the rriT spectra of pions and protons. We assume that pions and 
nucleons are in thermal equilibrium throughout the hydrodynamic evolution, and 
they freeze-out at a common temperature Tp — 120 MeV. We recall that there is no 
unique description so far for the pr spectra of hadrons. In [26] it has been shown 
that the 7r° spectra of WA80 and WA98 can be reproduced by perturbative QCD 
calculation if the initial px broadening is taken into account. Wang [26] has argued 
that the high pr 7r° spectra in central Pb -|- Pb collisions can not be due to collective 
flow. The pt broadening in heavy ion collisions can also be described reasonably 
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Figure 7: Same as fig. 6 for initial velocity profile of Eq. 6 with vo — 0.2. 
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Figure 8: The px distribution of pions for Pb + Pb collisions at CERN SPS energies. 
Left (Right): hadronic (QGP) initial state. The hadronic masses vary according to 
Eq. 15 with A = 1/2. The initial velocity profile is taken from Eq. 5 with vq = 0. 
The dotted curve in the left panel shows the pion spectra evaluated for the initial 
radial velocity of Eq. 5 with vq = 0.2. 



well by a random walk model [27] where transverse flow effects are not required 
to reproduce the data. Recently it has been proposed that the pr-broadening in 
high energy nuclear collisions can be generated by the Color Glass Condensate 
i.e. by the initial partonic phase formed after the collisions [28]. Keeping these 
possible scenarios of pr broadening in mind we do not attempt to reproduce the 
absolute normalization of the high px part of the pion spectra by hydrodynamic 
flow. We rather concentrate on the effects of the EOS containing the in-medium 
mass modification of hadrons on the pt spectra of pions and protons. We treat 
the normalization of pions and nucleons as parameters to be determined by the 
experimentally measured spectra. Negative hadrons and positive minus the negative 
are treated as pions and protons respectively. In fig. 8 the NA49 pion spectra is 
compared with both QGP and hadronic initial states. The hadrons (tt" and K~) 
from the resonance decays (p and 0) are about 5% of the direct pions at pr ~ 500 
MeV and are therefore neglected here. Results for QGP (hadronic) initial state is 
displayed in the left (right) panel in fig. 8. Both the initial states describe the data 
reasonably well for low transverse momentum. In fig. 9 the similar results are shown 
for protons. The experimental spectra is well reproduced by the hadronic as well as 
QGP initial state. However, due to various uncertainties which are discussed above 
it is not possible to state which of the EOS is realized in such collisions. 
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Figure 9: The pt distribution of protons for Pb + Pb collisions at CERN SPS 
energies. Left (Right): hadronic (QGP) initial state. The initial velocity profile is 
taken from Eq. 5 with vq — 0. The hadronic masses vary according to Eq. 15 with 
A = 1/2. 

In fig. 10 the Pt spectra for negative hadrons at RHIC energies are shown for 
the same value of the freeze-out temperature ( 120 MeV) and the EOS as shown 
in fig. 1. The agreement of our calculation with the experimental data [29] is 
reasonably well. In this case a thermalized QGP with initial temperature ~ 300 
MeV and thermalization time 0.5 fm/c is assumed. 

4 Summary and discussions 

We have solved the (3+1) dimensional hydrodynamical equations for QGP as well 
as hadronic initial states. The shift of hadronic masses at non-zero temperature are 
incorporated through the EOS. The FS is seen to be modified due to medium effects 
depending upon the magnitude of the shift in the hadronic masses. For a hadronic 
initial state in the universal scaling scenario (with the value of the exponent, A = 
1/2) the FS is larger compared to the case when no medium effects is considered. The 
difference in FS for QHD and vacuum scenario is neghgible due to small mass shift in 
QHD near the freeze-out point. The pion and proton spectra is well reproduced by a 
hadronic initial state when the mass shift of hadrons are taken into account via the 
EOS. It is seen that the pion and proton spectra is well reproduced with the same 
equation of state, initial conditions and the frcczc-out parameters which describe 
the WA98 [9] and the CERES/NA45 dilepton data [30], indicating a thermal source 
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Figure 10: The pt distribution of the pions at RHIC energies. The initial velocity 
profile is taken from Eq. 5 with = 0. The theoretical results are obtained for 
QGP initial state at a temperature 300 MeV. 

of initial temperature ~ 200 MeV at CERN SPS energies. 
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